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Active galaxies and gamma ray bursts are the sources of the highest energy photons detected 
by astronomical telescopes. We speculate that they may be the sources of the highest energy 
cosmic rays. This makes them true proton accelerators, where the highest energy photons 
are the decay products of neutral pions photoproduced when the proton beams interacts with 
ambient radiation. Neutrinos from the decay of charged pions represent an incontrovertible 
signature for proton acceleration. A main theme of this talk is that their fluxes can be 
estimated from the measured gamma ray luminosity by model-independent methods, based 
on dimensional analysis and textbook particle physics. 



Introduction 

Cosmic rays form an integral part of our galaxy. Their energy density is quali- 
tatively similar to that of photons, electrons and magnetic fields. It is believed 
that most were born in supernova blast waves. Their energy spectrum can 
be understood, up to perhaps 1000 TeV, in terms of acceleration by super- 
nova shocks exploding into the interstellar medium of our galaxy. Although 
the slope of the cosmic ray spectrum abruptly increases at this energy, par- 
ticles with energies in excess of 10^ TeV have been observed and, cannot be 
accounted for by this mechanism. The break in the spectrum, usually referred 
to as the "knee" can be best exhibited by plotting the flux multiplied by an 
energy dependent power E'^-'^^; see Fig. 1. 

The failure of supernovae to accelerate cosmic rays above 1000 TeV energy 
can be essentially understood on the basis of dimensional analysis. It is sensible 
to assume that, in order to accelerate a proton to energy E, the size R of the 
accelerator must be larger than the gyroradius of the particle in the accelerating 
field B: 
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Figure 1: Flux of high energy cosmic rays after multipUcation by a factor _E^-^^. Arrows 
point at structure in the spectrum near 1 PeV, the "knee" and 10 EeV, the "ankle'.' 



Therefore particles moving at the speed of hght c reach energies up to a max- 
imum value -Emax which must be less than 10"' TeV for the typical values of 
B and R for a supernova shock. Realistic modelling introduces inefficiencies 
in the acceleration process and, yields a maximum energy which is typically 
two orders of magnitude smaller than the value obtained by dimensional anal- 
ysis. One therefore identifies the "knee" with the sharp cutoff associated with 
particles accelerated by supernovae. 

Cosmic rays with energy in excess of lO^*' eV have been observed, some 
five orders of magnitude in energy above the supernova cutoff. Where cosmic 
rays with energy in excess of 1000 TeV are accelerated undoubtedly represents 
the most urgent problem in cosmic ray astrophysics and, one of the oldest un- 
resolved puzzles in astronomy. In order to beat the dimensional argument just 
presented, one has to accelerate particles over larger distances ^^^or identify 
higher magnetic fields B. This leads to a variety of speculations. EJ The large 
_B-fields in pulsars, near black holes and in the cores of active galactic nuclei, 
where they extend over large distances, have been exploited in models where 
the highest energy cosmic rays originate in point sources. 

Although imaginative arguments actually do exist to avoid this conclusion, 
it is generally believed that our galaxy is too small and its magnetic field 
too weak to accelerate the highest energy cosmic rays. Those with energy 
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in excess of 10^^ eV have gyroradii larger than our galaxy and should point 
back at their sources. Their arrival directions fail to show any correlation to 
the galactic plane, suggesting extra-galactic origin. As far as extra-galactic 
sources is concerned, nearby active galactic nuclei distant by ^ 100 Mpc are 
obvious candidates. The jets of blazars support magnetic fields of 10 G over 
distances of 10~^ parsecs or more. Using Eq. (3) we reach energies of 10^° eV, 
possibly higher because of beaming. Such speculations are reinforced by the 
fact that a cursory glance at the EGRET and Whipple results is sufficient to 
convince oneself that blazara^are also the dominant (exclusive?) sources of the 
highest energy gamma rays.a 

More "mundane" objects such as neutron stars and solar mass black holes 
may accelerate particle to 10^" eV in the most extreme circumstances: when 
a pair of them merges to produce a gamma ray burst. We return to this 
possibility at the end of this lecture. 



1 Active Galaxies as Proton Accelerators 

AGN are the brightest sources in the Universe. Their engines must not only 
be powerful, but extremely compact because their high energy luminosities are 
obserVjfid to flare by over an order of magnitude over time periods as short as 
a day.H Only sites in the vicinity of black holes, a billion times more massive 
than our sun, can possibly satisfy the constraints of the problem. Highly rel- 
ativistic and confined jets of particles are a common feature of these objects. 
It is anticipated that beams, accelerated near the black hole, are dumped on 
the radiation in the galaxy which consists of mostly thermal photons with 
densities of order lO^^/cm'^; see Fig. 2. The multi-wavelength spectrum, from 
radio waves to TeV gamma rays, is produced in the interactions of the acceler- 
ated particles with the magnetic fields and ambient light in the galaxy. In the 
more conventional electron models, the highest energy photons arc produced 
by Compton scattering of accelerated electrons on thermal UV photons which 
are scattered from 10 eV up to TeV energy D The energetic gamma rays will 
subsequently lose energy by electron pair production in photon-photon interac- 
tions with the radiation field of the jet or the galactic disk. An electromagnetic 
cascade is thus initiated which, via pair production on the magnetic field and 
photon-photon interactions, determines the emerging gamma-ray spectrum at 
lower energies. The lower energy photons, observed by conventional astro- 
nomical techniques, cirG, cis ci result of the cascade process, several generations 
removed from the primary high energy beams. 
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Figure 2: Possible blueprint for the production of high energy photons and neutrinos near 
the super-massive black hole powering an AGN. Particles, accelerated in sheet like bunches 
moving along the jet, interact with photons radiated by the accretion disk or produced by 
the interaction of the accelerated particles with the magnetic field of the jet. 



The EGRET instrument on the Compton Gamma Ray Observatory has 
detected high energy gamma-ray emission, in the range 20 MeV-30 GeV, from 
over 100 sources-B Of these sources 16 have been tentatively, and 42 sohdly 
identified with radio counterparts. All belong to the "blazar" subclass, mostly 
Flat Spectrum Radio Quasars, while the rest are BL-Lac objects.B In a unified 
scheme of AGN, they correspond to Radio_Loud AGN viewed from a position 
illuminated by the cone of a relativistic jet.l3 Moreover of the five TeV gamma- 
ray emitters identified by the air Chfirenkov technique, three are extra-galactic 
and are also nearby BL-Lac objectsJj The data therefore strongly suggests that 
the highest energy photons originate in jets beamed to the observer. Several of 
the sources observed by EGRET have^hown strong variability, by a factor of 2 
or so over a time scale of several days.EI Time variability is more spectacular at 
higher energies. On May 7, 1996 the Whipple telescope observed an increase 
of the TeV-emission from the blazar Markarian 421 by a factor 2 in 1 hour 
reaching, eventually, a value 50 times larger than the steady flux. At this point 
the telescope registered 6 times more photons from the MarkariaiJ-blazar, more 
distant by a factor 10^, than from the Crab supernova remnant.El 
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Does pion photoproduction by accelerated protons play a central role in 
blazar jets? This question has been extensively debated in recent yearsB If 
protons are accelerated along with electrons, they will acquire higher ener- 
gies, reaching PeV-EeV energy because of reduced energy losses. High energy 
photons result from proton-induced photoproduction of neutral pions on the 
ubiquitous UV thermal background. Accelerated protons thus initiate a cas- 
cade which dictates the features of the spectrum at lower energy .E£l From a 
theorist's point of view the proton blazar has attractive features. Protons, 
unlike electrons, efficiently transfer energy from the black hole in the presence 
of the high magnetic fields required to explain the confinement of the jets.EiJ 
Protons provide a "natural" mechanism for energy transfer from the central 
engine over distances as large as 1 parsec, as well as for the observed heating 
of the dusty disk over distances of several hundred parsecs.ll2l More to the 
point, the issue of proton acceleration can be settled experimentally because 
the proton bia^ar is a source of high energy protons and neutrinos, not just 
gamma rays.t2l 

Weakly interacting neutrinos can, unlike high energy gamma-rays and high 
energy cosmic rays, reach us from more distant and much more powerful AGN. 
It is likely that absorption effects explain why Markarian 421, the closest blazar 
on the EGRET list at a distance of 150 Mpc , produces the most prominent 
TeV signal. Although the closest, it is one of the weakest; the reason that 
it is detected whereas other, more distant, but more powerful, AGN are not, 
must be that the TeV gamma rays suffer absorption in intergalactic space 
by interaction with background infra-red lightJi3 This most likely provides 
the explanation why much more powerful quasars with significant high energy 
components such as 3C279 at a redshift of 0.54 have not been identified as TeV 
sources. Undoubtedly, part of the TeV flux is also absorbed on the infrared 
light in the source; we will return to this further on. 

First order Fermi acceleration offers a very attractive model for acceler- 
ation in jets, providing, on average, the right power and spectral shape. A 
cosmic accelerator in which the dominant mechanism is first order diffusive 
shock acceleration, will indeed produce a spectrum 

dN/dE (X E-"^ , 

with 7 ^ 2 -|- e, where e is a small number. For strong ultra-relativistic shocks 
it can be negative —0.3 ). Confronted with the challenge of explaining a 
relatively flat multi-wavelength photon emission spectrum which extends to 
TeV energy, models have converged on the blazar blueprint shown in Fig. 2. 
Particles are accelerated by Fermi shocks in bunches of matter travelling along 
the jet with a bulk Lorentz factor of order 7 ~ 10. Ultra-relativistic beaming 
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with this Lorcntz factor provides the natural interpretation of the observed 
superluminal speeds of radio structures in the jet.Ej In order to accommodate 
bursts lasting a day in the observer's frame, the bunch size must be of order 
TcAt ^ 10^^ parsecs. Here F is the Doppler factor, which for observation 
angles close to the jet direction is of the same order as the Lorentz factor.LI 
These bunches are, in fact, more like sheets, thinner than the jet's width of 
roughly 1 parsec. The observed radiation at all wavelengths is produced by the 
interaction of the accelerated particles in the sheets with the ambient radiation 
in the AGN, which has a significant component concentrated in the so-called 
"UV-bump'.' 

In electron models the multi-wavelength spectrum consists of three com- 
ponents: synchrotron radiation produced by the electron beam on the i?-field 
in the jet, synchrotron photons Compton scattered to high energy by the elec- 
tron beam and, finally, UV photons Compton scattered bw the electron beam 
to produce the highest energy photons in the spectrum. □ The seed photon 
field can be either external, e.g. radiated off the accretion disk, or result from 
the synchrotron radiation of the electrons in the jet, so-called synchrotron-self- 
Compton models. The picture has a variety of problems. In order to reproduce 
the observed high energy luminosity, the accelerating bunches have to be po- 
sitioned very close to the black hole. The photon target density is otherwise 
insufficient for inverse Compton scattering to produce the observed fiux. This 
is a balancing act, because the same dense target will efficiently absorb the 
high energy photons by 77 collisions. The balance is difficult to arrange, espe- 
cially in light of observations showing that the high energy photon flux extendja 
beyond TeV energy .0 The natural cutoff occurs in the 10-100 GeV region.ta 
Finally, in order to prevent the electrons from losing too much energy before 
producing the high energy photons, the magnetic fleld in the jet has to be 
artificially adjusted to less than 10% of what is expected from equipartition 
with the radiation density. 

For these, and the more general reasons already mentioned in the introduc- 
tion, the proton blazar has been developed. In this model protons as well as 
electrons are accelerated. Because of reduced energy loss, protons can produce 
the high energy radiation further from the black hole. The more favorable 
production-absorption balance far from the black hole makes it relatively easy 
to extend the high energy photon spectrum above 10 TeV energy, even with 
bulk Lorentz factors that are significantly smaller than in the inverse Compton 
models. Two recent incarnations of the proton blazar illustrate th at th ese mod- 
els can also describe the multi-wavelength spectrum of the AGN.EjEfl Because 
the seed density of photons is still much higher than that of target protons, 
the high energy cascade is initiated by the photoproduction of neutral pions by 
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accelerated protons on ambient light via the A resonance, The protons collide 
either with synchrotron photoiia, produced by electrons,Ej or with the photons 
radiated off the accretion disk,E3 as shown in Fig. 2. 

2 The Neutrino Flux from Blazar Jets 

Model-independent evidence that AGN are indeed cosmic proton accelerators 
can be obtained by observing high energy neutrinos from the decay of charged 
pions, photoproduced on the A resonance along with the neutral ones. The 
expected neutrino flux can be estimated in six easy steps. 

1. The size of the accelerator R is determined by the duration, of order 1 
day, over which the high energy radiation is emitted: 

R = Ttc= 10"^ parsecs for t = 1 day. (4) 

2. The magnitude of the i?- field can be calculated from equipartition with 
the electrons, whose energy density is measured experimentally: 

- — = ^(electrons) ^ 1 erg/cm . (5) 
2/io 

This yields a value for the magnetic field of 5 Gauss. A similar value is 
obtained by scaling _B-fields in the jets of Fannaroff-Riley type II galaxies 
at kiloparsec distances, to the Markarian 421 luminosity, and to trans- 
verse distances in the milliparsec rangcO 

3. In shock acceleration the gain in energy occurs gradually as a particle 
near the shock scatters back and forth across the front gaining energy 
with each transit. The proton energy is limited by thCj-Iifetime of the 
accelerator and the maximum size of the emitting regior£3 R, 

E<KZeBRc. (6) 

Here Ze is the charge of the particle being accelerated and B the ambi- 
ent magnetic field. The upper limit basically follows from dimensional 
analysis. It can also be derived from the simple requirement that the 
gyroradius of the accelerated particles must be contained within the ac- 
celerating region R. The numerical constant K ^ Q.\ depends on the 
details of diffusion in the vicinity of the shock, which determine the ef- 
ficiency by which power in the shock is converted into acceleration of 
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particles. In some cases it can reach values close to 1. The maximum 
energy reached is 

i^max = eBRc = 5 X 10^^ eV 

for B = 5 Gauss and R = 0.02 parsecs. We here assumed that the boost 
of the energy in the observer's frame approximately compensates for the 
efficiency factor, i.e. KT ~ 1. 

The neutrino energy is lower by two factors which take into account i) 
the average momentum carried by the secondary pions relative to the 
parent proton {{xp) ~ 0.2) and ii) the average energy carried by the 
neutrino in the decay chain tt"*" i^fil^~^ c^v^v^, which is roughly 1/4 
of the pion energy because equal amounts of energy are carried by the 
four leptons. The maximum neutrino energy is 

-B.max = ^max {xp) \ ^ 10^^ eV , (7) 

i.e. neutrinos reach energies of 10^ PeV. 

4. The neutrino spectrum can now be calculated from the observed gamma 
ray luminosity. We recall that approximately equal amounts of energy 
are carried by the four leptons that result from the decay chain 7r+ — > 
t'^/i"'" e^v^Vfj^. In addition the cross sections for the processes p7 — > 
p7r° and P7 mr~^ at the A resonance are in the approximate ratio of 
2:1. Thus 3/4 of the energy lost to photoproduction ends up in the 
electromagnetic cascade and 1/4 goes to neutrinos, which corresponds to 
a ratio of neutrino to gamma luminosities {L^, : L^) of 1 : 3. This ratio is 
somewhat reduced when taking into account that some of the energy of 
the accelerated protons is lost to direct pair production (p+7 — > e+e~p): 

In order to convert above relation into a neutrino spectrum we have to 
fix the spectral index. We will assume that the target photon density 
spectrum is described by a power law, where a is small for AGN 

with flat spectra. The number of target photons above photoproduction 
threshold grows when the proton energy Ep is increased. If the protons 
are accelerated to a power law spectrum with spectral index 7 (= 2 + 
e), the threshold effect implies that the spectral index of the secondary 
neutrino flux is also a power law, but with an index flattened by (1 + a) 
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as a result of the increase in target photons at resonance when the proton 
energy is increased: 



dE„ 



E„ 



-(l+e-Q) 



(9) 



For a standard non relativistic shock with e = and a flat photon target 
with a = 0, the neutrino spectrum will flatten by just one unit giving 
E^^ ~ constant. From Eqs. (5) and (6) 
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The calculation is stable as long as e — a is smaller than 1 because the 
luminosity integral is not sensitive to the lower limit of the integration. 

5. Assuming that the high energy 7 ray flux from Markarian 421 results from 
cascading of the gamma ray luminosity produced by Fermi a|Ccelerated 
protons, we obtain the neutrino flux from the measured valuen of of 
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where the numerical estimate corresponds to a = e 



(11) 

and the value 



of Ey max of Eq. (4) . This calculation reveals that for the small values 
of e and a anticipated, the neutrino flux is essentially determined by the 
value for E^, max- 

In order to calculate the diffuse flux from the observed blazar distribu- 
tion, we note that the EGRET collaboration has constructed a luminosity 
function covering the observation of the ^^20 most energetic blazars and 
estimated the diffuse gamma ray luminosity .1111 From the ratio of the 
diffuse gamma ray flux and the flux of Markarian 421, we obtain that 
the effective number of blazars with Markarian 421 flux is ~ 130 sr~^. 
The diffuse neutrino flux is now simply estimated by multiplying the 
calculated flux for Markarian 421 by this factor. A correction for the 
difference in spectral indices of gamma ray and neutrino fluxes enhances 
the neutrino flux by a factor of three. The flux corresponds to an energy 
regime well below the high energy cut-off. The transition to the cutoff 
should be smooth because of the superposition of the different redshifts 
and cut-off energies of the individual blazars. 
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This concludes our calculation. It illustrates how the proton blazar, un- 
like the electron blazar, requires no large Doppler factors and no fine-tuning 
of parameters. For the proton blazar, radiation and magnetic fields are in 
equipartition, the maximum energy matches the BR value expected from di- 
mensional analysis and, finally, the size of the bunches is similar to the gyrora- 
dius of the highest energy protons. It is not a challenge to increase gamma ray 
energies well beyond the TeV energy range. Reasonable variations of the val- 
ues of magnetic field strength B, the efficiency parameter K and the Doppler 
boost factor T may allow us to account for the highest energy cosmic rays with 
E r^3x 10^° eV. 

Also, our calculation demonstrates why modclSS which differ in many 
aspects, yield very^imilar predictions for the neutrino flux, consistent with the 
ones obtained heral3; this is illustrated in Fig. 3. 
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Figure 3: Diffuse neutrino flux from blazars. The numerical result of Eq. (11) multipliaiiihTi 
3 X 130 sr~^ and corrected for redsliift in the cutoff is compared to recent calculations.lijil3 
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3 Cosmic Rays and High Energy Neutrinos 



Rather than scahng to the TeV gamma ray flux, we can use the cosmic ray flux 
at ultra high energies to braclcet expectations for the neutrino flux. Models for 
proton acceleration in hot spots of Fanaroff-Riley tmpe II galaxies can explain 
the observed cosmic ray spectrum above ~ 10^* eVJ13 This energy corresponds 
to the "ankle" in the spectrum, where the observed spectral index flattens 
from 3 to 2.7. The model requires an E~'^, or flatter, injection spectrum which 
steepens above 10"'^^ eV to the observed E^"^-"^ spectrum as a result of energy 
loss in the source- interactions with the microwave background, and cosmo- 
logical evolutionJi^ Because of the strict limitations on the density of target 
photons at the acceleration site, previpusly discussed, roughly similar neutrino 
and proton luminosities are expected.ES In order to understand this balance it 
is important to realize that in astrophysical beam dumps the accelerator and 
production target form a symbiotic system. Although larger target density 
may produce more neutrinos, it also decelerates the protons producing them, 
in a delicate acceleration-absorption balance. Equal cosmic ray and neutrino 
luminosity implies: 

j dE^{E^ dN^/dE^) - LcK ^ 10'^ TeV cm^^ sr^^ . (12) 

The bulk cosmic ray luminosity has been conservatively estimated by assuming 
that it is due to an E~^-'^ spectrum above ^ 10^^ eV. This spectrum has been 
normalized to the observed EeV cosmic rays. It is interesting to note that this 
luminosity is. a factor 25 below the measured diffuse gamma ray luminosity 
from AGN.tll This is, within an order of magnitude, in agreement with the 
relation of neutrino and gamma ray luminosities of Eq. (5), and, if anything, 
implies a conservative estimate of the neutrino flux. 

Assuming a generic E~'^ neutrino spectrum, the equality of cosmic-ray and 
neutrino luminosities implies: 

^ dN, 10-10 _2 -1 -1 no^ 
'''^dE^^EA^)'^ ' ■ (^3) 

A not too different result is obtained by assuming equal numbers of neutrinos 
and protons, rather than equal luminosities. It is clear that our estimate is 
conservative because the proton flux reaching Earth has not been corrected for 
absorption of protons in ambient matter in the source, or in the interstellar 
medium. 



11 



4 Event Rates in Underground Muon Neutrino Telescopes 

The probability to detect a TeV neutrino is roughly 10^^.0 It is easily com- 
puted from the requirement that, in order to be detected, the neutrino has 
to interact within a distance of the detector which is shorter than the range 
of the muon it has produced. In other words, in order for the neutrino to be 
detected, the produced muon has to reach the detector. Therefore, 



where is the muon range and Ajnt the neutrino interaction length. For 
energies below 1 TeV, where both the range and cross section depend linearly 
on energy, n — 2. Between TeV and PcV energies n = 0.8 and A = 10^^, with 
E in TeV units. For EeV energies n = 0.47, A = 10"^ with E in EcV. 

We are now ready to compute the diffuse neutrino event rate by folding 
the neutrino spectrum of Eq. (8) with the detection probability of Eq. (11). 
We also multiply by 130 sr^^ for the effective number of sources: 



which implies a yield of two neutrinos every three days in a kilometer-scale 
detector, assuming only 27r coverage. 

The steeper, but lower luminosity, flux of Eq. (10) predicts more events 
when folded with Eq. (11), about 150 km~^ year~^ sr~^ assuming that the flux 
extends down to TeV energy. The result does not depend strongly on the lower 
limit of the neutrino integral, it only drops by a factor of three if the neutrino 
flux flattens below 100 TeV_We again conclude that a kilometer-scale neutrino 
detector may be requiredEj It is however important to realize that, had we 
assumed a E~^ spectrum, the resulting flux would have scaled with the ratio 
of luminosities to about an order of magnitude below Eq. (12). The energy 
dependence of the detection efficiency of underground muon neutrino detectors 
is such that most of the events are detected in the high (low) energy end for a 
E^^ [E'"^) spectrum. 

5 Evidence for the Proton Blazar? 

Astronomy with protons becomes possible once their energy has reached a 
value where their gyroradius in the microgauss galactic field exceeds the di- 
mensions of the galaxy. Provided that intergalactic magnetic fields are not 
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too strong, protons with 10^° eV energy point at their sources with degree- 
accuracy. At this energy their mean-free-path in the cosmic microwave back- 
ground is unfortunately reduced to only tens of megaparsecs. A clear window 
of opportunity emerges: Are the directions of the cosmic rays with energy in 
excess of ~5 x 10^^ eV correlated to the nearest AGN (red-shift z less than 
0.02), which are known to be clustered in the so-called "super-galactic" plane? 
Although far from conclusive, there is some evidence that such a correlation 
may exist.EJ Lack of statistics at the highest energies is a major problem. 
Future large aperture cosmic ray detectors such as the new Utah HIRES air 
fluorescence detector and the Auger giant air shower array will soon remedy 
this aspect of the problem.!!! 

Another problem is that the pointing accuracy is not really understood. It 
depends on the distance d to the source and the gyroradius in the intergalactic 
magnetic field: 

d dB 



i?gyro E 



or scaled to units relevant to the problem 

9 ^ ( 1 Mpc) (io-3g) 



0.1° / E 



10^" cV 



(16) 
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Speculations on the strength of the inter-galactic magnetic field range from 
10^^ to 10^^^ Gauss. For a distance of 100 Mpc, the resolution may therefore 
be anywhere from sub-degree to nonexistent. It is reasonable to expect that 
magnetic fields are higher in regions where matter is clustered. Higher val- 
ues may therefore be appropriate for the local cluster and the super-galactic 
plane, regions relevant to this problem. Optimistically, one may anticipate 
that future high statistics experiment such as HIRES in Utah and the Auger 
giant air shower array will provide indirect information on the magnitude of 
the magnetic fields between galaxies. 

We have already drawn attention to the 10 TeV maximum photon en- 
ergy as the demarkation line between the electron and proton blazars. The 
~ 10 GeV cutoff in the inverse Compton model can be pushed to the TeV 
range in order to accommodate the Whipple data on Markarian 421, but not 
beyond. Bringing the accelerator closer to the black hole may yield photons in 
excess of 10 TeV energy — they have, however, no chance of escaping without 
energy loss on the dense infrared background at the acceleration site. HEGRA 
has been monitoring the 10 closest blazars, including Markarian 421, with its 
dual telescope systems: the scintillator and the naked photomultiplier detec- 
tor arrays. The announcementJ23 that their upper limit on the photon flux of 
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50 TeV and above for the aggregate emission from the ten nearest blazars, 
may be a signal, could provide the first compelling evidence that blazar jets 
are indeed proton accelerators. 

In summary, there are hints that active galaxies may be true particle ac- 
celerators with proton beams dictating the features of the spectrum. With the 
rapidly expanding BaikaLand AMANDA detectors producing their first hints 
of neutrino candidates j2jO observation of neutrinos from AGN would establish 
the production of pions and identify the acceleration of protons as the origin of 
the highest energy photons. A definite answer may not be known until these 
detectors reach kilometer size. Neutrino telescope builders should take note 
that, although smaller neutrino fluxes are predicted than in the generic AGN 
models of a few years ago,Ea they are all near PeV energy where the detection 
efficiency is increased and the atmospheric neutrino background negligible. Be- 
cause of the beaming of the jets, the neutrinos have a flatter spectrum peaking 
near the 10^ TeV maximum energy. The actual event rates are, in the end, 
not very different. 

If confirmed, these models strongly favor the construction of neutrino tele- 
scopes following a distributed architecture, with large spacings of the optical 
modules and relatively high thresholdBj This also opens up opportunities for 
alternative techniques such as the radio technique, or the detection of hor- 
izontal air showers with giant air shower arrays.cj Optimists, on the other 
hand, can find reasons to anticipate the discovery of AGN neutrinos with much 
smaller telescopes. With a sufficiently high proton target density in the accel- 
eration region, much larger fluxes of neutrinos may be produced in a^proton- 
proton cascade. The predicted fluxes are however model-dependent £3 It is 
also possible, even likely, that accelerated protons which produce neutrinos do 
not escape the source, or escape after significant energy loss. Such absorption 
effects increase the neutrino flux relative to the observed high energy cosmic 
ray flux, also leading to larger neutrino fluxes. 

6 Comments on the Highest Energy Cosmic Rays: Particles with 
Macroscopic Energy 

Ever since the pioneering Haverah Park experiment0 revealed that cosmic 
particles are accelerated up to 10^° eV energy, the origin of these very highest 
energy particles has been hotly debated. The recent observation of isolated 
events with even higher energy is nothing less than paradoxical; they imply 
aspects of particle physics or astrophysics not revealed in previous experiments. 

In October 1991, the_Ely's Eye cosmic ray detector recorded an event of 
energy 3.0±q;54 x 10^'' eV.ca This event, together with an event recorded by the 
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Yakutsk air shower array in May 1989^ of estimated energy ~ 2_>< 10^" eV, are 
the two highest energy cosmic rays ever seen. A recent paperuJ reports that 
the Akeno Giant Air Shower Array, an instrument of over 100 scintiUation 
detectors spread over a 100 km^ area in the Kofu vineyards in Japan, also 
recorded an event of (1.7-2.6) x 10^" eV energy. These cosmic rays' energies 
exceed by a factor of a hundred miUion those of the highest energy particle 
beams achieved with man-made accelerators. When colliding with atmospheric 
nuclei, the center of mass energy is of the order of 700 TeV or ~ 50 Joules, 
almost 50 times that of the future Large Hadron Collider at CERN. 

Nature accelerates microscopic particles to macroscopic energy by mecha- 
nisms that are still a matter of speculation, e already noted that the gyrora- 
dius of a 10^^ eV proton in the 3 x 10"^ gauss galactic field is roughly 10,000 
light-years, comparable to the size of our galaxy. So, 10^° eV particles should 
travel in straight paths from their sources through the galactic and intergalactic 
magnetic fields. Textbook particle physics is sufficient to derive that 10^" eV 
protons travelling through the omnipresent 3K photon background will pho- 
toproduce pions, and will thus be demoted in energy over a distance of less 
than 10 Mpc, i.e. much less than the 100 Mpc plus distance from the posited 
sources. Although one can do astronomy with the highest energy charged par- 
ticles, their reach into the Universe is very limited by attenuation in the cosmic 
photon background. 

Here are the facts. The mean free path of a 3 x 10^" eV proton in the 
cosmic photon background is only 8.8 Mpc. The probability for a proton to 
traverse 100 Mpc without an interaction is 1.16 x 10~^. Alternatively, a cosmic 
ray proton needs an energy of 3 x 10^^ eV to reach Earth from a 100 Mpc source 
with the observed energy. Needless to say, acceleration to IQ^^ eV becomes 
a challenge, even in the context of supermassive black holesrd In the case of 
the Fly's Eye event no obvious candidate sources exist near or within the error 
box of the observed arrival direction. (There are some candidate AGN sources 
for which the distance is, however, unknown). From the previous discussion it 
is clear that the identification of the highest energy cosmic rays with protons 
is problematic. The above arguments apply, mutatis mutandis, to nuclei. 

As previously mentioned, above a threshold of about 5 x 10^^ eV, the 
so-called Greisen-Zatsepin cutoff, cosmic rays lose energy in the universal pho- 
ton background by photoproducing pions which eventually decay into photons 
and neutrinos. It has therefore been anticipated that at the highest energies, 
7-rays and neutrinos dominate the spectrum. We already discussed that ac- 
tive galaxies may be powerful neutrino sources. Weakly interacting neutrinos 
escape acceleration sites without absorption and travel without attenuation 
through interstellar space. They may form the dominant high energy compo- 
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neiit of cosmic ray accelerators if the highest energy protons (and photons) 
are absorbed in the ambient matter in the source or in intersteUar spa|Cfi_The 
proton origin of the high energy events has therefore been questioned. E3o 

Photons of 10^° eV energy interact with the Earth's magnetic field by pair- 
producing electrons and start a cascade well befoiie entering the atmosphere. 
The highest energy cosmic rays are not photonsH Neutrino origin is, on the 
contrary, consistent with the observed air shower profile. A very high energy 
muon-neutrino interacting high in the atmosphere transfers roughly one fifth of 
its energy to an air nucleus which subsequently produces a hadron-like cascade. 
Such an event would, of course, be atypical since the atmosphere is relatively 
transparent to the weakly interacting neutrinos. It would be the tip of an 
iceberg that reaches its maximum inside the Earth, waiting to be discovered 
by the high energy neutrino telescopes now under construction. 

One may argue that proton and neutrino assignments are equally prob- 
lematic, by a factor 10~^. The factor represents the probability of reaching 
Earth through the microwave background in the case of a proton and the ra- 
tio of the neutrino to proton interaction cross section with the atmosphere in 
the case of neutrino origin. So, the skeptic will conclude that these particles 
are neither protons, nor photons nor neutrinos, and that nothing accelerates 
particles to 10^° eV anyway. This is the paradox. The optimistic theorist may 
argue that we are observing a new particle, created in a topological defect, ef- 
ficiently decelerated from Planck energies. The experimentalist will argue that 
we need more data. This is quite a challenge given a flux of order 10"'^/ (km^ 
sr year), which points at the need for a detector with an effective area of 
order 10000 km^ which will be hopefully commissioned as the "Auger Giant 
Air Shower Detector'.' In the meantime Utah's HIRES has a shot at the same 
problem. 

7 Are GRBs the Sources of the Highest Energy Cosmic Rays? 

Our speculation that AGN jets are the sources of the highest energy cosmic 
rays has been very much inspired by the fact that they are also the sources 
of the highest energy photons. Photons of GeV energy are also produced in 
gamma ray bursts (GRBs)Efl and this raises the possibility that two of the 
outstanding enigmas in astronomy, high energy cosmic rays and GRBs, may 
actually represent a single puzzle. 

GRBs, typically, release a fraction of a solar mass of energy 10^'^ ergs) 
over millisecond time scales in photons with a very hard spectrum. Approach- 
ing supernova-scale energy release over such short time scales, GRBs are clearly 
not in the realm of the stars and plasmas of everyday astronomy. They must 
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represent true high energy cosmic accelerators, possibly merging neutron stars 
or black holes. 

As we have seen, similar speculations emerge when faced with the formidable 
challenge p£ identifying the origin of the highest energy cosmic rays. It has been 
suggestedE3 that, though unknown, the same cataclysmic mechanism that pro- 
duces GRBs also produces the highest energy cosmic rays. This association 
is reinforced by other features of the data besides the phenomenal energy and 
luminosity of their respective sources: 

• with improved statistics it has become increasingly suggestive that both 
GRBs and the highest energy cosmic rays are produced in cosmological 
sources, i.e. distributed throughout the Universe, and 

• the average rate (over volume and time) at which energy is injected 
into the Universe as gamma rays from GRBs is similar to the rate at 
which energy is injected in >10^'^ eV cosmic rays in order to produce the 
observed cosmic ray flux. 

The obvious question is whether GRBs and cosmic rays are actually ob- 
served in coincidence? As noted above, the high energy protons lose energy on 
the universal background of cosmic photons, the highest energy cosmic rays 
observed at Earth reach us from no further than a few tens of megaparsecs. 
It is an observational fact that the rate of GRBs within this limited cosmo- 
logical volume is only one every few hundred years, while 10^° eV cosmic rays 
are detected at a rate of one every few years. There is no contradiction here! 
As a result of deflection and energy dispersion of the cosmic rays by random 
magnetic field the source pulse is time broadened by non-straight, lengthened 
paths. Intergalactic magnetic fields of order 10~^^ Gauss are sufficient to 
spread the cosmic ray signal over 100 years. This, unfortunately, implies that 
one does not expect a sharp correlation between the observed directions and 
arrival times unless GRB are not truly cosmological and their concentration is 
somehow "locally" enhanced within a radius of tens of megaparsecs. 

Independently of the magnitude of the magnetic fields in the source and 
in the interstellar medium, sharp coincidences between photons and charged 
cosmic rays are excluded because of propagation effects inside our own galaxy. 
The deflection angle of a 10^" eV cosmic ray is expected to be of order 1 degree 
as a result of propagation over kiloparsec distances in the microgauss magnetic 
field of our galaxy. The angle is inversely proportional to the energy of the 
particle. The wriggling of the trajectory of the cosmic ray by the magnetic field 
produces delay in its arrival of order several months, scaling inversely with the 
square of the energy. As previously discussed, cosmic ray times are expected 
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to have been spread by hundreds of years rather than months. This has not 

discouraged the search for possible coincidences. The Utah Fly's Eye detector 
observed the highest energy cosmic ray ever on October 15, 1991. Within the 
error box of its direction the BATSE detector on the Compton Gamma Ray 
telescope observed GRB (910503) 5.5 months earlier. This event happened to 
be the highest fluency event in the BATSE catalogue. This burst also had 
photons with energy as high as 10 GeV. With the usual pitfalls of post-factum 
statistics it is difficult to judge the weight of evidence of such association. Is 
this an atypical nearby event? Are the GRBs not truly cosmological? Is this 
association fortuitous? The future will tell. Other, less dramatic associations 
of GRBs and high energy cosmic rays have been suggested. Independently of 
the polemic, cosmic rayers will undoubtedly look out for future events in the 
direction of the February 17, 1994 "Superbowl" event with a fluency twice that 
of GRB (910503). 

Particularly intriguing is that the highest energy particles are emitted with 
delays of tens of minutes. It is possible that the delayed photons observed in 
events like GRB 940217, where a 25 GeV photon arrived with a delay of 77 min- 
utes, are the decay products of neutral pions photoproduced by accelerated 
protons interacting with the MeV photons in the GRB shock. This assump- 
tion is sufficient to compute the number of high energy neutrinos produced in 
the shock. 

Prom the observed (and theoretically expected) high energy GRB gamma 
ray spectrum 

we can fix the photon flux by normalizing to 1 event above 25 GeV, per EGRET 
area of 0.1 m^ and per burst. We obtain 

250 

= — GeV m-2 burst"^ . (19) 

Pion photoproduction is dominated by the A resonance. Particle physics re- 
lates the neutrino and gamma ray luminosities by Eq. (8). This is a conserva- 
tive estimate because it assumes no high energy gamma ray absorption in the 
source. Eqs. (8) and (19) determine the neutrino flux. 

The probability to detect a neutrino can now be computed from (14) and 
(15). For energies below 10'^ GeV, where both the range and cross section 
depend linearly on energy, n = 2 and A ~ 10"^^. Above this energy we simply 
assume n = 0, which is adequate for this calculation. The number of detected 
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events is obtained by folding the neutrino flux with the detection probability 



(20) 



Here E^^th is the threshold of the detector. For the flux of Eq. (19) (multiplied 
by 3/13) we obtain 
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dE 

E^ 



Therefore 



5 X 10"* m"^ burst" 



10 km ^ year 



(21) 



(22) 



For 1 burst per day one predicts more than 10 events per squared kilometer per 
year. These events are relatively easy to detect because they come at precise 
times from known directions, supplied by the gamma ray observations. One 
may detect that many events from a single burst that is 10 times closer, and, 
possibly, a handful of events in the operating AMANDA detector with an area 
of ~ 10'' m^. Notice that the muons are near TeV energy provided the detector 
threshold is below this value, as is the case for AMANDA with a threshold of 
~ 50 GeV. 

The calculation can be done in a slightly different way, yielding the same 
result. We can think of the GRB phenomenon as a source of energy releasing 
a sizeable fraction of a solar mass every day. This corresponds to a flux of 
10^'' ^ 10^* ergs/sec. For an E~'^ spectrum Eq. (20) can be conveniently 
rewritten in the form 



= 10 



]^q38 ej-g/s 



10 kpc 



(lO"' m^ year) 



(23) 



i.e., 10 neutrinos of TeV energy are observed in an AMANDA-size detector per 
year for a galactic source with Eddington luminosity. For the GRB sources 
the flux is increased l)y ahnost ten orders of magnitude, though the distance 
d ~ 3 X lO"' Mpc {{z) ^ 1) is more than 10'"' times larger. The increased 
luminosity does not quite compensate for distance to the sources and Eq. (23) 
yields 10 ~ 10^ events per year in a kilometer squared detector. This conflrms 
our previous estimate. 

In the previous estimates we have obviously assumed that kinetic energy in 
the shock is efflciently converted into the acceleration of protons, an assump- 
tion intrinsic to models where GRBs are the sources of the highest energy 
cosmic rays. Our calculations assume that pions are the source of the observed 
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gamma xays. While such modela^J may be somewhat disfavored, Waxman and 
BahcallEZI have pointed out that the production of high energy neutrinos may 
be a feature of more conventional models as well. In their model the gamma 
rays are produced by synchrotron radiation by relativistic electrons accelerated 
in the shock, possibly followed by inverse-Compton scattering. Provided that 
the efficiency with which kinetic energy is converted to accelerated protons 
is comparable to that for electrons, the GRBs can still generate the highest 
energy cosmic rays.Ej These protons will photoproduce pions and neutrinos on 
the MeV gamma rays in the burst. The predicted rates are similar to the ones 
estimated above, but the neutrinos are some two orders of magnitude higher 
in energy. The predicted neutrino flux is 

dN/dE = A/E^ for E > E^ (24) 

= A/{Ei,E) for E < Eb (25) 

with A = 4 X 10^13 XeV (cm^ s si)-^ and E^ ~ 100 TeV. 

The flux is similar to the one obtained for AGN in the previous section. 
This is not surprising as the mechanisms are similar. It is higher than the AGN 
flux, but most neutrinos are produced near the maximum energy of 100 TeV, 
rather than 100 PeV. The event rates are in the end similar. 

There is also a very different class of models where the delayed high energy 
gamma rajp^are produced in the interactions of fireball protons with hydrogen 
gas clouds. Ea For the case of photoproduction, pions arc photoproduced on the 
A resonance and textbook particle physics implies that the neutrino luminosity 
is 3/13 of the gamma ray luminosity. For this alternative model where the 
gamma rays originate in nuclear interactions, the luminosities of photons and 
neutrinos should be roughly equal, and, if anything, the neutrino luminosity 
larger. The question is whether absorption effect (of gammas in the source) 
are small enough for these predictions to survive. 

The association of GRBs and high energy cosmic rays has other powerful 
observational implications. Models predict fluxes of high energy photons and 
neutrinos which exceed cosmic ray proton fluxes by a couple orders of magni- 
tude and I guess that, in the near future, air Cherenkov telescopes, air shower 
arrays, neutrino telescopes and, of course, the 10000 km^ "Auger" cosmic ray 
detector will begin to probe this intriguing association. 
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